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ABSTRACT 

Chandra observations of low-luminosity supermassive black holes in nearby elliptical galaxies provide 
tight limits on both their nuclear luminosities and on their Bondi accretion rates. We examine Chandra 
constraints on NGC 6166, the dominant galaxy in the cluster Abell 2199, which hosts a ~ 10 9 M Q 
black hole. We measure a nuclear X-ray luminosity Lx.ikcV ~ 10 40 ergs _1 and show that the density 
and temperature profiles of the hot interstellar medium imply a Bondi accretion rate of M Bondi ^ 
3 x 1CU 2 Mq yr _1 . This accretion rate predicts a nuclear luminosity of ~ 10 44 ergs _1 for a canonical 
radiative efficiency of 10%. Unless the Bondi estimate is inappropriate and/or the accretion rate onto 
the black hole is significantly reduced, the observed nuclear flux constrains the radiative efficiency of the 
accretion flow to be rj ss 10 -5 . We show that low-radiative efficiency accretion flows with radial density 
profiles p °S r~ 3 / 2 (and not significantly flatter) can explain the observed nuclear X-ray luminosity but 
that the power output from the jets in NGC 6166 is also important to the energetics of the system. 



1. INTRODUCTION 

Recent dynamical measurements of black hole masses in 
nearby galaxies have transformed black hole studies, allow- 
ing us to quantify directly their expected accretion power. 
In those cases in which the fuel supply to the black holes 
can be estimated, it has become apparent that the ex- 
pected power far exceeds the observed luminosity. This 
situation is perhaps best illustrated by the case of nearby 
elliptical galaxy nuclei (e.g. Fabian & Rees 1995; Di Mat- 
teo et al. 2000a, hereafter DMOOa and references therein) 
and our Galactic center (e.g. Narayan et al. 1998 and 
references therein). There are only two possible explana- 
tions for the low-luminosities of nearby black holes: (a) 
the accretion occurs at extremely low rates or (b) the ac- 
cretion occurs at low radiative efficiencies as predicted, for 
example, by advection dominated accretion models (e.g. 
Rees et al. 1982; Narayan & Yi 1994,1995; Abramowicz et 
al. 1995). 

Assuming that the accretion in elliptical galaxies occurs 
primarily from the hot, quasi-spherical interstellar medium 
(ISM), Bondi theory (1952) can be used to estimate the 
accretion rates onto the supermassive black holes. Such es- 
timates require accurate measurements of both the density 
and temperature of the hot ISM at the 'Bondi accretion 
radius', the radius where the gravitational potential of the 
central black holes begins to dominate the dynamics of the 
hot gas. 

In order to unambiguously determine whether the low 
luminosities of nearby black holes are due to a low radia- 
tive efficiency in the accreting gas, it is also necessary to 
measure the nuclear power. When combined with the esti- 
mated accretion rates, this gives us a direct measurement 
of the radiative efficiency of accretion, rj. 

Thanks to its high spatial resolution and sensitivity, the 
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Chandra X-ray observatory is able, for the first time, to 
detect nuclear X-ray point sources in nearby galaxies and 
provide us with direct measurements of their luminosities. 
Chandra also allows us to measure the central densities 
and temperatures of the ISM close to the accretion radii of 
the central black holes, and therefore determine the Bondi 
accretion rates in these systems to much greater accuracy 
than before. 

In this Letter we explore the implications of Chandra 
observations of the giant elliptical galaxy NGC 6166. We 
report the detection of a nuclear point source in the galaxy 
with a luminosity L x ~ a few 10 40 ergs _1 . We measure 
the central X-ray gas density and temperature and calcu- 
late the Bondi accretion rate. We show that if the central 
black hole of ~ 10 9 M is fed at the estimated Bondi rate, 
the inferred accretion radiative efficiency is rj < 10~ 5 . We 
show that at the given accretion rate ADAF models can 
explain the observed nuclear luminosity. The presence of 
outflows in the accretion flows is also consistent with the 
present constraints. We discuss possible reasons why the 
Bondi accretion rate estimate might be inappropriate and 
show that a significant contribution to the nuclear flux is 
likely to be due to emission from the jets. 



2. NUCLEAR LUMINOSITY AND ACCRETION RATE 

NGC 6166 (z = 0.031) is a cD galaxy at the center of 
the cooling flow cluster A2199, and is therefore located in a 
highly gas-rich environment. The central black hole mass 
is estimated by Ferrarese & Merritt (2000) using the cor- 
relation between the masses of black holes and the velocity 
dispersions of their host bulges (Mbh oc ct 4 7 ). For NGC 
6166, er c = 326 km s -1 which gives Mbh ~ 10 9 M Q , im- 
plying a discrepancy of about an order of magnitude with 
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the value obtained by Magorrian et al. (1998), who esti- 
mated M BH ~ 10 10 M Q . 

NGC 6166 is known to host an active nucleus, classified 
as an FR I source, which powers two symmetric parsec- 
scale radio jets and radio lobes. The jets are close to 
the plane of the sky (Giovannini et al. 1998). Beaming 
should not therefore be important. Given the measured 
black hole mass, the Eddington luminosity of NGC 6166 
is ^Edd ~ 10 47 ergs _1 . However, the observed bolometric 
luminosity is < 10 41 ergs _1 . 

We adopt the cosmological parameters Hq = 
65kms~ 1 Mpc~ , = 1 and A = 0, which imply a dis- 
tance to NGC 6166 of 142 Mpc). 



The Chandra observations of the A2199 cluster were car- 
ried out on 1999 December 11 using the Advanced CCD 
Imaging Spectrometer and back-illuminated CCD detec- 
tors. The net exposure time was 14.95 ks. Here, we only 
highlight results for the central, dominant galaxy NGC 
6166; details of the cluster and point source analysis are 
discussed by Johnstone et al. (2000). 

The Chandra data for NGC 6166 reveal the presence of a 
nuclear point-like X-ray source (detected over background 
at 7 sigma significance) with a surface brightness distribu- 
tion consistent with the point-spread function (PSF) of the 
instrument (a ~ 0.45 arcsec). The X-ray source is coinci- 
dent with the peak of the emission from the central radio 
source (Giovannini et al. 1998). In order to accurately de- 
termine the central point source flux, we have modeled the 
central surface brightness profile using a simple power-law 
(7x oc r~ a ; accounting for the underlying ISM) plus Gaus- 
sian model, convolved with the instrumental PSF. We find 
that the surface brightness profile within the innermost 10 
kpc radius region can be well-described by a power-law 
model of slope a = 0.4 ± 0.025 and a point source with 
55 ± 27 counts (la errors; where the error in the counts is 
dominated by the uncertainty in the slope of the density 
profile) . 

We have extracted a spectrum for the central source us- 
ing a region of ~ 1.5 arcsec radius (which should contain 
> 90 per cent of the flux from the source) and used a 
background spectrum taken from the 1.5 — 2.5 arcsec an- 
nulus. Given the faintness of the central source and the 
moderate exposure time of the observation, the spectral 
parameters are not firmly constrained and the data can 
be described by either a power law or a thermal plasma 
model (mckal). The power-law model gives a photon in- 
dex T = 1.54+ 0^5 assuming Galactic absorption, and a 
flux density at 1 keVof 7 ± 2 x 10~ 15 ergcm~ 2 s" 1 keV -1 , 
where all parameters are quoted at a la confidence. (Al- 
lowing for intrinsic absorption, the maximum intrinsic col- 
umn density is A^ ~ 2 x 10 21 atomcm~ 2 for T = 2.9 
and F(lkeV) = 14 x 10~ 15 ergenr 2 s" 1 keV -1 ; see Fig. 3 
and contour plots in Johnstone et al. 2000). The flux de- 
termined from the spectrum is consistent with the value 
measured from the surface brightness profile, assuming a 
consistent spectral model for the nuclear emission. 

2.2. M 

We can estimate the accretion rate from the hot ISM 
onto the central black hole using Bondi accretion the- 



ory. Matter passing within the accretion radius R A ~ 
(c/c s ) 2 Rs is assumed to be accreted. Here c is the speed of 
light, c s ~ 10 4 T 1 / 2 cm s _1 is the sound speed (with T the 
gas temperature), and i?s = 2GM/c 2 the Schwarzschild 
radius of the black hole. 

We have measured the temperature and density profile 
of the intracluster gas within the central 100 kpc radius 
region. Direct measurements of these quantities are only 
possible down to within r ~ 1 arcsec of the nucleus, due 
to contamination from the nuclear point source. The filled 
triangles in Fig. 1 show the temperature profile determined 
by fitting the annular counts with a single temperature 
thermal emission model (mekal) acted on by Galactic and 
(free-fitting) intrinsic absorption. The open squares show 
the deprojected temperature profile in which the inner re- 
gions arc corrected for the emission from further out seen 
in projection. The temperature is higher in the outer re- 
gions, where the cluster potential dominates and decreases 
(and flattens off) in the inner regions where the galaxy po- 
tential takes over. 

The density profile shown in Fig. 2 was derived by de- 
projecting the surface brightness profile (e.g. Fabian et al. 
1981) and matching the deprojected temperature profile. 
The density profile (Fig. 2) is steeper in the outer regions 
(r > 10 kpc) and follows the canonical cooling-flow profile 
of p 25 r~ x . Inside this radius the density profile flattens to 
a value p oc r -°- 45 ± - 09 ( we exclude the central 1.5 arcsec 
region, which is affected by the point source). 

We take kT ~ 1.3 keV as measured in the innermost 
regions. This gives 

R A ~ 0MT{-}M 9 kpc, (1) 

where Mg = 10 9 M Q which corresponds to to ~ 0.08 arc- 
sec. 

The accretion rate is related to the density and temper- 
ature by the continuity equation, with the velocity v = c s ; 

M = 4irR 2 A p A c s (R A ) (2) 

where pa and c s (Ra) are the density and sound speed at 
the accretion radius. 

Given that the density and temperature cannot be mea- 
sured directly at the accretion radius we assume kT = 
1. 3+0^7 keV and estimate the density at Ra by extrapolat- 
ing the best-fit power-law model (see above) which gives 
n(R A ) ~ 0.6 ±0.1 cm" 3 . 

This implies 

M Bo „di = 1.31J;? x 10 24 Ml n . 6 gs" 1 (3) 

- 3 x 10~ 2 Meyr 1 , 

(where we have neglected the uncertainties on Mbh)- 
Eqn. (2) gives an estimate of the accretion at the outer 
edge of the accretion flow. As discussed in the next ses- 
sion, the mass accretion rate onto the black hole may be 
smaller if e.g. it decreases with radius because of an out- 
flow. 

This Bondi accretion rate implies L Bondi = 
V^BondiC 2 ~ 10 44 ergs _1 if t] = 0.1, as in a standard 
radiatively efficient thin disk. This is about four orders 
of magnitude greater than the observed luminosity of the 
central X-ray source, implying that the radiative efficiency 
has to be rj ~ 10~ 5 suggesting that hot accretion flows 
with low radiative efficiencies may be relevant. 
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3. MODELS 

In a hot accretion flow around a supermassive black hole, 
the majority of the emission arises in the radio and X-ray 
bands. In the radio band the emission results from syn- 
chrotron radiation from the inner parts of the accretion 
flow. The X-ray emission is due to either bremsstrahlung 
emission or inverse Compton scattering of the soft syn- 
chrotron photons (e.g. Narayan, Barret & McClintock 
1998) 

We measure radii in the flow in Schwarzschild units: 
R = rR s , where R s = 2GM BU /c 2 . We measure black 
hole masses in solar units and accretion rates in Edding- 
ton units: Mbh = mM Q and M — mMEdd- We take 
M Edd = 10L E dd/c 2 = 2.2 x KT 8 mM Q yr" 1 , i.e., with a 
canonical 10% efficiency. We take r = 10 4 to be the outer 
radius of the flow. The Bondi accretion rate in Eddington 
units is meondi = 1-3 x 1CP 3 . 

The predicted spectrum from an ADAF depends 
(weakly) on the ratio of gas to magnetic pressure (3, the 
viscosity parameter a, and the fraction of the turbulent 
energy in the plasma which heats the electrons, S. Here, 
we fix a = 0.1, (3 = 10, and 5 = 0.1. The two major pa- 
rameters, though, are the accretion rate M and the black 
hole mass Mbh, both of which are constrained. With 
Mbh assumed, we normalize the models to the observed 
Chandra flux. This gives us the m required by the models 
to explain the X-ray emission. This has to be consistent 
with our estimate of meondi- 

The solid line in Fig. 3 shows the predicted spectrum 
for a pure inflow ADAF model fitted to the 1 keV flux. 
In this model m ~ 10~ 3 , consistent with the Bondi esti- 
mate. Comptonization of the synchrotron emission dom- 
inates the X-ray emission in this model. ADAFs mod- 
els can therefore explain the observed X-ray emission. At 
the given Bondi rate the radiative efficiency of an ADAF 
is ~ 10~ 5 corresponding to that derived in §2.2. Note, 
however, that in these models the synchrotron component 
cannot explain the VLBI radio emission (Giovannini et al. 
1998; Fig. 3). This implies that the majority of the con- 
tribution in the radio band is likely to be due to emission 
from the jets (e.g., dotted line). 

Following the proposal of Blandford & Begclman (1999) 
that mass loss may be important in hot accretion flows, 
Di Matteo et al. (1999, hereafter DM99; DMOOa) and 
Quataert & Narayan (1999; QN99) have constructed mod- 
els in which a significant fraction of the accreted mass in 
an ADAF is lost to an outflow rather than being accreted 
onto the central object. In that work (see also Di Matteo, 
Fabian & Carilli 2000b, hereafter DMOOb) it was shown 
that such models are favored (or required) over the pure 
inflow ADAF for explaining the observed spectral energy 
distributions of elliptical galaxy nuclei. 

The importance of outflows in ADAFs can be parame- 
terized by using different radial density profiles in the flows 
with p cx r ~ 3 / 2 +P and < p < 1 (i.e. flatter than in the 
pure Bondi inflow with p oc r~ 3 / 2 ). This implies mass in- 
flow rates which satisfy the relation M = M out (r/r out ) p , 
where we should satisfy M ou t ~ Meondi and r out = 10 4 
as before. Bremsstrahlung emission in these models dom- 

2 Note that in a CDAF the Bondi accretion rate estimate is not ap 
rate at which matter is fed into the flow. The accretion rates (at large 



inated the X-ray emission. This produces a harder X-ray 
spectrum, most of which (in the 0.5 - 10 keV band) is 
produced in the outer regions of the flow (the tempera- 
tures in the inner regions approach 100 keV to 1 MeV) at 
f ~ ^out where m ou t ~ m (see details in DMOOa; QN99). 
The dashed line shows a model with p = 0.1 which implies 
m ou t = 3 x 10~ 3 . If p > 0.1, TO out becomes inconsistent 
with the Bondi estimate (Eq.3) unless the errors on Mbh 
are sufficiently large (e.g., for p — 0.2, m ou t ~6x 10~ 3 , so 
Mbh would need to have been underestimated by a factor 
> 50%). 

For certain regimes (e.g. for a < 0.1) convection might 
also become important in hot accretion flows, (e.g. Stone, 
Pringle & Begelman 1999; Quataert & Gruzinov 2000; 
Narayan, Igumenshchev & Abramowicz 2000). Although 
CDAFs have a very different structure than ADAFs, they 
have density profiles p cx r -1 / 2 (corresponding to the 
p = 1 case of the outflow models above). Ball, Narayan 
& Quataert (2000) have constructed spectral models for 
CDAFs similar to those considered by DM99 and QN99 
for outflows 2 . Taking p = 1, even with m ou t ~ meondi, 
implies very small contributions from a CDAF to the X- 
ray flux. A CDAF if present cannot explain the observed 
emission. 

Finally, the dotted line in Fig. 3 shows an example of 
non-thermal jet models for the full band spectrum. It is 
clear that the observed radio jets must contribute signif- 
icantly to the observed radio and optical emission even 
if an ADAF (with or without outflows) is also present. 
Fig. 3 shows that the full spectral energy distribution 
of NGC 6166 can be easily explained by the emission 
from different compact regions emitting self-absorbed non- 
thermal synchrotron radiation. The typical parameters 
for these components are B ~ 7 x 10~ 4 ^ghz^ 1 bi 2 G and 
t t ~ 0.4T| 12 i/ tG H Z ; where T B12 = 10 12 K is the bright- 
ness temperature and v t the self-absorption frequency (see 
Di Matteo, Carilli & Fabian 2000 for details, hereafter 
DMOOb). Emission from a jet is far more efficient even 
if it originates from relatively low density regions and can 
easily give rise to the observed X-ray flux (see also e.g. Fal- 
cke 1999). Explaining most of the observed emission as 
being due to the jets would be consistent with the unify- 
ing model in which FR I radio galaxies are mis-oriented 
BL Lac objects and the fact that the radio-to-optical and 
optical-to- X-ray spectral indexes of NGC 6166 are con- 
sistent with those of radio selected BL Lac objects (see 
Fig. 3 in Hardcastle & Worrall 1999). However, this in- 
terpretation would still not relax the requirement for a 
low-radiative efficiency in the accretion flow. 

4. DISCUSSION 

We have examined Chandra constraints on the low 
luminosity active nucleus in NGC 6166. By adopting 
Mbh — 10 9 M Q and inferring the ISM density and tem- 
perature at i?A we have estimated the Bondi accretion 
rate to be ~ 3 x 10~ 2 M Q yr _1 . This value is < 10 times 
lower than the typical Bondi rates determined in previous 
studies based on ROSAT High Resolution Imager (HRI) 
observations of other nearby elliptical nuclei (e.g., M87, 
NGC4649, NGC4472; DMOOa) and is more consistent with 

propriate as convection and possibly conduction significantly alter the 
radii) arc expected to be lower than in the ADAFs 
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the upper limits inferred from high frequency radio stud- 
ies of the accretion rates at small radii in the ellipticals 
(DM99;DM00b; Wrobel & Herrnstein 2000). 

The Chandra observations of NGC 6166 indicate a nu- 
clear luminosity of L x _ . 5-7 koV ~ 4xl0 40 ergs _1 . Thisim- 
plies that if matter is accreted from the ISM at the Bondi 
rate the radiative efficiency of the accretion flow has to be 
r] < 10 -5 . We have shown that this is consistent with the 
predictions of pure inflow ADAF models. Hot accretion 
flow models with significantly flatter density profiles (with 
strong outflows or CDAFs), which are favored on theoret- 
ical grounds and from radio observations (DM99,DM00b), 
cannot explain the observed nuclear flux (unless Mbh is 
significantly larger). This may imply that most of the ob- 
served flux is dominated by non-thermal emission from the 
jets. 

Although with Chandra observations we can infer Bondi 
accretion rates with much greater confidence than before, 
for NGC 6166 we can only measure the ISM densities 
and temperatures directly at a radius ~ 1 kpc which is 
still a factor ~ 15 larger than Ra (§2.2). Given the dis- 
tance to NGC 6166, previous ROSAT HRI observations 
would have only been able to measure the density pro- 
file down to r > 10 kpc, where we observe the density 
profile to steepen. Extrapolation of this steeper density 
profile to the accretion radius would then have lead to 
higher inferred accretion rates. Even with Chandra, how- 
ever, we can still only infer the ISM density at Ra in NGC 
6166, which implies that we could still be overestimating 
the Bondi accretion rate (if e.g. the density profile flat- 
tened further towards Ra)- Given that in order to explain 



the observed nuclear flux, an ADAF requires accretion 
rates similar to those currently estimated, a decrease in 
the Bondi estimate would make the presence of an ADAF 
undetectable and their emission negligible. Similar argu- 
ments apply if the temperature closer to Ra increased. 

It is also possible that the Bondi calculation itself may 
be inappropriate. If strong convection were important, as 
in CDAFs, the rates at which matter is fed through the 
accretion radius could be significantly reduced. Heating of 
the ISM by the jets (DMOOb) could also decrease signif- 
icantly the Bondi accretion radius and hence the rate at 
which mass is fed into the accretion flows. If the accretion 
rate is low enough, the requirement for low radiative effi- 
ciency in the accretion flow may be significantly relaxed. 

We have shown that emission from the jets in these 
low-luminosity nuclei provides an important contribution 
to the overall energy output of NGC 6166. Determin- 
ing the variability properties of the X-ray emission may 
be important for resolving the jet contribution from that 
of an ADAF. Chandra observations of other, closer, low- 
luminosity black holes (e.g., M87) are needed to further 
assess the question of low radiative-efficiency versus low- 
M accretion in nearby black holes. 

We thank Ramesh Narayan for many useful comments. 
T. D. M. acknowledges support for this work provided by 
NASA through Chandra Postdoctoral Fellowship grant 
number PF8-10005 awarded by the Chandra Science Cen- 
ter, which is operated by the Smithsonian Astrophysical 
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Fig. 1. — The temperature profile in NGC 6166 measured by Chandra. The open squares show the deprojected temperature 
measurements and the filled triangles the observed (projected) values. 
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Fig. 2. — The deprojected density profile in the inner regions of A 2199. 
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Fig. 3. — Spectral models calculated for hot accretion flows fitted to the Chandra X-ray flux. The solid line is an ADAF 
model. The required accretion rate, tin = 8 x 10~ 4 , is consistent with the inferred Bondi value. The dashed line shows a 
model with outflows. Here p = 0.1, m out = 3 x 10~ 3 . The dotted line is the sum of self-absorbed, non-thermal synchrotron 
emitting regions in the jets. Beaming is ignored since the jets lie very close to the plane of the Sky (Giovannini et al. 1998). 
The solid dots are VLBI flux measurements from Giovannini et al. (1998), the open circles the VLBI core peak fluxes and 
the HST (Chiaberghe et al. 1999) continuum flux and ROSAT HRI upper limit (Hardcastle & Worrall 1999). The bow-tie 
shows the Chandra flux (in the 0.5 — 7 kcV band) for the range of power-law slopes with the hatched area corresponding 
to the ones acceptable with Galactic absorption fixed. 



